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Abstract —We study the performance of centralized least mean- 
squares (CLMS) algorithms in wireless sensor networks where 
nodes transmit their data over fading channels to a central 
processing unit (e.g., fusion center or cluster head), for parameter 
estimation. Wireless channel impairments, including fading and 
path loss, distort the transmitted data, cause link failure and 
degrade the performance of the adaptive solutions. To address 
this problem, we propose a novel CLMS algorithm that uses a 
refined version of the transmitted data and benefits from a link 
failure alarm strategy to discard severely distorted data. Fur¬ 
thermore, to remove the bias due to communication noise from 
the estimate, we introduce a bias-elimination scheme that also 
leads to a lower steady-state mean-square error. Our theoretical 
findings are supported by numerical simulation results. 

Index Terms —centralized parameter estimation, LMS adaptive 
algorithms, wireless sensor networks, fading channels. 

I. Introduction 

ENTRALIZED least mean-squares (CLMS) algorithms 
were introduced in previous studies for applications in 
networks with multiple processing nodes JT|, J2)- In the 
development of these algorithms, it is commonly assumed that 
the communication links between the nodes and the fusion 
center are error free, i.e., there is no channel distortion such 
as fading, path loss, and noise during the data transmission. 
Consequently, these algorithms may experience substantial 
performance degradation when used over wireless sensor net¬ 
works (WSN), as evidenced in recent studies on distributed 
estimation in WSN operating over fading channels |3l, Pfl. 

Reference a and El propose novel distributed adaptive 
algorithms for networks with imperfect communication links. 
However, these algorithms do not include provisions to handle 
possible link failures as well as the estimation bias resulting 
from the link noise. Work in |5) studies the bias compensation 
problem for scenarios where only the regression data are noisy. 
Previous work (6l investigates the problem of space-varying 
parameter estimation over networks using distributed process¬ 
ing approaches. In these works, however, the communication 
links between the nodes are assumed to be ideal, i.e., no fading 
and communication noise are taken into consideration. 

We, in this paper, bring all these factors into account and 
study the performance of CLMS algorithms for parameter 
estimation in WSN in the presence of radio channel im¬ 
pairments. Specifically, we propose a novel CLMS algorithm 
that uses a refined version of the transmitted information and 
benefits from a link failure alarm strategy to discard severely 
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distorted data. We demonstrate the presence of a bias in the 
adaptive parameter estimates, consequence of the link noise, 
and introduce a bias-elimination scheme that also significantly 
decreases the steady-state mean-square error (MSE) of the 
network. Beside the application of our proposed bias-removal 
technique in centralized networks, it can be also used in 
networks with a distributed hierarchical structure, where the 
cluster heads act as fusion centers and distributed processing 
is performed at a cluster level. 

The adaptive feature of the proposed algorithm is highly 
desirable in non-stationary signal environments where the 
underlying network parameters change over time. In summary, 
the contribution of this letter are as follow: a) development 
of a new CLMS algorithm for parameter estimation in WSN 
operating over fading channels; b) analysis of the developed 
algorithm and identification of the main technical issue under 
such condition, namely the estimation bias; c) derivation of the 
relationship between the bias estimates with that of the channel 
perturbations; d) development of a new bias-compensation 
technique to remove the bias and e) performance analysis 
of the developed bias-compensated algorithm, including the 
derivation of the stability conditions and the steady-state mean- 
square deviation (MSD) expressions. 

II. Problem Formulation 
Consider a network of N randomly distributed sensor nodes 
that are deployed over a geographical area to estimate an 
unknown parameter vector w° £ C Mxl . At time instant 
i€N, node k £ {1,2, - - - ,7V} collects a measurement datum, 
d k {i) £ C, that is related to w° through a linear modeQ. i.e.: 

dk(i ) = Uk,i,w° + v k {i) (1) 

where the regression vectors u k ,i £ C lxM are zero-mean 
i.i.d. in time and independent over space with positive def¬ 
inite covariance matrix R Utk = E[iz£ iU k ,i\, the noise terms 
v k {i) £ C are zero-mean i.i.d. in time and independent over 
space with variance k . The regression data u kA and the 
measurement noise v/(j) are independent (T), [2), 13- 
Let us assume that each node k assembles its available data 
into a row vector as z k ,i = [u k ,i, d k (i)] and sends it over 
a wireless channel to the network fusion center. During its 
transmission, the data experience distortion of the form: 

Zk,i = a k (i)z kti + (2) 

where z k ,i = [uk,i, d k (i)\ denotes the received data at the 
fusion center and v' k \ £ C lx ( M + 1 ) j s the link noise with 

*In this letter, the same mathematical notation as in 0 is used, where 
boldface fonts are reserved for random variables. 
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v kl = l v k l 'l v k d \ i )\- We consider a k (i) = h k {i)J ^ for 
analog data transmission, where h k {i) £ C represents the 
fading channel coefficient, P is the transmit signal power, r k 
is the distance between node k and the fusion center and a is 
the path loss exponent. For digital-type data communication, 
the coefficient a k {i) will be an indicator function that models 
the link-failure (packet loss) in baseband after decoding. 
The fading coefficients h/Ji) are zero-mean circular complex 
Gaussian, i.i.d. in timqj and independent over space with 
variance a 2 k . The link noise vectors v^ k \ are zero-mean, 
i.i.d. in time and independent over space with covariance 
matrix u^ 2 /, while the components v k l j and are 

independent. The random variables h k (i), u k ,i, v k (i ) and v ki 
are all mutually independent. We use r° k to denote the free- 
space transmission range of node k, as obtained using Friis 
formula for a the given power budget, antenna gains and carrier 
frequency 0. Furthermore, we let <; k represent the threshold 
signal-to-noise ratio (SNR) of the received signal from node 
k , defined as the SNR of the received signal over a non-fading 
link with communication range r k , i.e., q k = P/cr[^ k (r k ) a . 
We use <; k (i) to denote the link instantaneous SNR, which is 
given by <; k (i ) = \h k {i)\ 2 P/cr^ k r k . The transmission from 
node k to the fusion center will be successful if <;/.('/) exceeds 
the threshold level q k , i.e., if \h k (i)\ 2 > [r k /r k ) a \ otherwise, 
the link fails. In the digital case, this amounts to a packet loss, 
while in the analog case this would correspond to a temporary 
loss of synchronization. Since the channel coefficients, h k (i), 
are zero mean circular complex Gaussian random variables, 
their squared magnitudes, \h k (i)\ 2 , have exponential distribu¬ 
tion ID. Let A = l/a 2 k , then the probability of successful 
transmission will be: 

Pk = Pr(|M*)| 2 > (pr) a ) = e- x ^f/ r ^ a (3) 


III. Centralized LMS (CLMS) Algorithm 

Let us first consider the analog-type data communication 
between nodes and the fusion center. In this case, the fu¬ 
sion center pre-process {u k ,i, dk{i)} k =i to partially recover 
{u k p, d k (i)} k=1 before using them in the LMS iteration. This 
pre-processing can be of linear form, e.g.: 

u k ,i = g k {i)u k ,i (4) 

d k (i) = g k {i)d k {i) (5) 


where g k (i) are scalar equalization coefficients. Assuming 
negligible channel estimation error, these coefficients can be 
obtained, in practice by, e.g., the least squares (LS) method: 


9 k (i) 


K{i) [rj[ 
IM*)| 2 V P 


( 6 ) 


Substituting <(6l) into ® and (0 leads to: 

Uk,i = u k ,i+g k (i)v ( k J (7) 

d k [i) = d k {i) +g k (i)v ( k \i) (8) 


The network can now seek the unknown parameter vector w° 
using the pre-processed data by minimizing: 


lr rhis assumption is true when the time interval between successive itera¬ 
tions of the adaptive process is larger than the coherence time of the channels. 


N 


k= 1 


J(w) = E Ik (*) I dk (i) - u k ,iW | 


(9) 


where X k {i) is a random variable with Bernouilli distribution 
defined as: 


Z k (i) 


1, if Sk(i) > 

0, otherwise 


( 10 ) 


That is, Xfe(z) = 1 when the transmission from node k to the 
fusion center is successful and Ik(i) = 0 otherwise. We note 
that the probability of success, i.e., p k = E[Xfc(z)], is given by 
0. Since the cost function @ is strictly convex, its optimal 
point will be its only stationary point. This leads to ED: 

N _i N 

It^ctrl = f ^ ^ (p k Ru,k T ^ ^ Pfc T d.Utk') (11) 

fc= 1 fe=1 


where r dUtk = E [d k (i)u* k4 \, and 

R Vtk =E[l k (i)\g k (i)\ 2 ]R%l (12) 

with R^fl = E[v k ‘f*v k 1 ^]. The expectation term on the right 
hand side of (IT2l) can be obtained as: 

rVk r a . I r-g a \ 

E[X k (i)\g k (t)\ 2 ] = / -^Xe X ^~ Xk >dy k (13) 

where y k = ^ p and x k = 

Since the covariance and cross-covariance of the data may 
not be available in practice, the estimate ic" r | given by (ITTI) can 
alternatively be sought using the following CLMS algorithm: 

N 

Wi = Wi -1 + /J,^z k (i)u* k ^d k (i) - u k ,iWi- 1 ) ( 14 ) 

k= 1 

where p > 0 is the step-size. The CLMS algorithm for digital 
data transmission model, i.e., when a k (i ) = X k (i) will be 
the same except that u k ,i and d k {i) are now replaced by 
Uk,i and d k (i), respectively. In this case, X k (i) will be an 
indicator function that models the link-failure or packet loss 
due to fading. We now proceed to investigate the convergence 
of the CLMS algorithm. To simplify the presentation, we first 
consider the case a k (i) = 1 and assume that the links fail with 
probability p k . Later, we include fading and path loss effects 
into the analysis. Under this condition, from 0, we have: 

d k (i) = d k (i) + v ( k \i) (15) 

Uk,i = Uk,i + V^j (16) 

Using 0. (00 and (fl6l >. we then obtain: 

d k (i) = u k ,iW° + v k {i) (17) 

Vk{i) = v k (i) + v ( k \i) - v ( k }w° (18) 

The variance of v k {i) is then given by = a 2 k + a^ k + 

w 0 *R^lw°. We now define the weight error-vector Wi = w°— 
Wi and use (fl4] > and ( fTTI) to arrive at: 

w i = D i Wi^ 1 -pt i (19) 


where D, = I - and — 

^2 k= i X k [i)u* k iV k (i). Taking expectation from both sides of 
<03. we obtain: 

E[zhi] = UEfzhj-i] — ptf (20) 
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where D = E[D., ; ] = I - pJ^k=iPk(Ru,k + R^ k ) and t f = 

E[U] = — X^fcLi PkR^fl w °- From recursion (120b , it can be 
verified that the CLMS algorithm will be mean-stable if 
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0 < /r < 


Amax^X)fc=l Pk(Ru,k + 


( 21 ) 


If the step-size satisfy OH . from (l20ll . we can also show that 
the mean estimate of the algorithm deviates from the optimal 
estimate by: n 

b = w°- _lim E[wi] = ~2pk(R u ,k + R^lj) t f ( 22 ) 
k= 1 


will be the same if we replace u kl and d k (i) with u kl and 
d k (i), respectively and set g k (i) to one. 

We now proceed to analyze the proposed algorithm. Similar 
to previous section, we first consider the case with a k (i ) = 1. 
Under this condition, it can be verified that the weight error 
vector of the algorithm evolves with time according to: 

Wi = DiWi -1 - gti (26) 

where, in this case, £)* = I-gJ2k=i (*)(«£,»«*:,* ~R^l) 

and U = Y^k=i X k{i)(u* k4 v k (i) + R^l™ 0 )- By taking the 
expectation of (f26l> . we arrive at: 


We observe that the bias, b, does not depend on the step- 
size /r. Therefore, reducing the step-size will not improve the 
accuracy of the estimated parameters. 

To incorporate the effects of fading and path-loss into the 
above results, we can use expressions © and ® and repeat 
the analysis. Doing so, we obtain similar results as before, 
except that in this case the power of the communication noises 
and increase by the factor of E[Xfc(z)|g fe (i)| 2 ]. 

Consequently, <EU and (Ell can still be used if, in these 
expressions, we replace PkR^l with R Vl k- 

IV. Bias-Compensated CLMS 

As seen from ( 1221 . the mean of the weight error vector of 
the CLMS algorithm in ( 1 1 41 . i.e., E[wr,j, converges to a non¬ 
zero vector. In what follows, we propose a bias-compensation 
scheme and develop a new form of centralized LMS algorithm 
whose mean estimation error converges to zero. 

As shown in d22li . the estimation bias b is due to the term 
tf, which in turn is caused by the regression noise component 
v k J j. Assuming zero regression noise in ©, the unbiased 
optimal estimate of the network will be: 

n « n 

w° = (j2PkRu,k) (5 ~^Pkfdu,k ) (23) 

fc=1 fc=1 


E[«tj] = BE[ti)j_i] (27) 

where D = E[£),] = I — P^Z k=l PkRu,k- It can be verified 
that lim^oo E[iL;] —>• 0 if 

0 < fi < -t —y -y (28) 

^max ( —i PkRu.k ) 

The mean-square deviation (MSD) recursion of this algorithm 
can be expressed as: 

E||«ti||| = ElltUi-ilIfy + (vec(/2 2 X)) *vec(£) (29) 

where E\\ibi\\% = E||w° - Wi\\%, E > 0, S' = E [D.^D.,] 

and 

N 

x = J2p k (<k f Ru,k + (P - 1) (Ri:l^°w 0 *R^l) 

k =1 

+ R^lTr(w°w°*R^l)) (30) 

Using this recursion, the steady-state MSD of the algorithm 
converges to: 

7ctri = [vec {ji 2 X)]*(I - F)- 1 vec(/ M ) (31) 

where F ~ F) T V) D* . As mentioned before, the steady-state 
MSD of the network under fading condition can be computed 
by replacing PhR^'i with R Vtk given by (|T2l ). 


We, in this paper, propose a CLMS algorithm that runs over 
the noisy received data and achieves the optimal estimate (1231 . 
The basic idea in our development is to construct an objective 
function whose gradient vector is identical to that of the cost 
© with zero regression noise iTffl . The following objective 
function satisfies this criterion: 

N 

J(w) = £)E[X fc (i)(|d fc (*) - u k ,iW | 2 - \\R^ k w \| 2 )] (24) 

fc=l 

It can be verified that the Hessian of J ( w ) is positive definite 
and hence this cost function is strongly convex. From (l24l >. 
we then arrive at the following bias-compensated CLMS (BC- 
CLMS) algorithm where p, is the new step-size: 

N 

Wi = + p'^X k {i)^u* ki (d k {i) - Uk,iWi- 1 ) 

k= 1 

+ 9k(i)R%l w i- 1 } ( 25 ) 

This algorithm requires the covariance matrix of the regression 
noise, which can be estimated off-line or during the operation 
of the algorithm in real-time as in e.g., 0, G2- The bias- 
compensated CLMS algorithm for digital data transmission 


V. Simulation Results 

In this section, we present computer experiments to illustrate 
the performance of the CLMS algorithm (fl4l i and the BC- 
CLMS algorithm d25l > over a wireless sensor network. We con¬ 
sider a WSN consisting of N = 5 nodes, uniformly distributed 
over a squared area of [x , y] = [1,1] km with the fusion center 
located at x = y = 0.5km. The unknown parameter vector to 
be estimated is w° = [—0.8006, —0.3203 +j0. 1601,0.4804] T 
with M = 3. All nodes have equal transmit signal power P = 
lOmW, the nominal transmission range is r k = 0.3km, and the 
path loss exponent is a = 2.5. We initialize the algorithms 
with tu_i = [0,0,0] T and set the step-size to p = 0.003. 
To generate communication noise vectors v^ k \, we adopt 
complex normal distributions with zero mean and variances 
cr^l 2 e {0.0617,0.0560,0.0923,0.0831,0.0476}. The mea¬ 
surement noise v k (i) is generated using a normal distribution 
with variances ct 2 k € {0.069,0.090,0.087,0.092,0.061,}. 
The fading coefficients of the radio links are obtained from 
a circular complex normal distribution with zero mean and 
unit variance. The regression data u k l are generated from 
a circular complex Gaussian distribution where the (n,m)- th 
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Fig. 1: Mean error vector versus time for ( |25V 
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Fig. 2: MSD versus time for CLMS (14) and BC-CLMS (25). 
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Fig. 3: The MSD performance of BC-CLMS algorithm under 
channel estimation error and temporal correlation. 


remains invariant. In contrast, as the channel estimation error 
increases the steady-state MSD of the algorithm increases. As 
this result indicates, the proposed algorithm is robust to small 
channel estimation error and temporal correlation. 


VI. Conclusion 

We studied the performance of CLMS algorithms for pa¬ 
rameter estimation over wireless sensor networks where links 
between nodes and the fusion center are impaired by fading 
and noise. The analysis and the numerical simulations show 
that the proposed BC-CLMS algorithm is stable and converges 
to an unbiased estimate for small adaptation step-sizes. 


entry of correlation matrix R u is given by for m > n , 

with = k/(2N). Figure Q] shows the mean convergence 
performance of the BC-CLMS algorithm (1251) . In this figure, 
we illustrate the convergence of the real part of the mean 
error vector As the results show the mean error 

vector converges to zero after 1200 iterations. In contrast, 
the mean error vector of the CLMS algorithm (fl4l > converges 
to the biased value of b = —0.1489 — j'0.0086, —0.0278 + 
jO.0317, 0.1064 — j0.0025] T , which is in agreement with our 
analytical finding. We examine the mean-square performance 
of CLMS and BC-CLMS algorithms for two cases. In the 
first case, the fading coefficients change at each time i and 
links fail according to the channels instantaneous SNR. In 
our results, we refer to this case as the fading case. For 
the second case, the communication links are distorted by 
additive noise but links fail with probability of 1 — pk with 
Pk = [0.342,0.336,0.415,0.363,0.473]. We refer to the sec¬ 
ond case as the link-failure case. Figure [2] shows the MSD 
of the algorithms for these two cases. As expected, the BC- 
CLMS algorithm achieves a significantly lower steady-state 
MSD than the CLMS in both cases. In Fig. [3] we compare the 
performance of the proposed BC-CLMS algorithm with that 
of the conventional CLMS algorithm proposed in J2)- As this 
figure shows the proposed algorithm outperforms CLMS with 
about 8dB. This figure also shows the performance of the BC- 
CLMS algorithm under imperfect channel information at the 
fusion center and under temporal channel correlation. We ob¬ 
serve that as channel correlation increases the converges speed 
of the algorithm decreases while its steady-state performance 


References 

[1] C. G. Lopes and A. H. Sayed, “Diffusion least-mean squares over 
adaptive networks: Formulation and performance analysis,” IEEE Trans, 
on Signal Processing, vol. 56, no. 7, pp. 3122-3136, July 2008. 

[2] F. S. Cattivelli and A. H. Sayed, “Diffusion LMS strategies for dis¬ 
tributed estimation,” IEEE Trans, on Signal Processing, vol. 58, no. 3, 
pp. 1035-1048, Mar. 2010. 

[3] R. Abdolee and B. Champagne, “Diffusion LMS algorithms for sensor 
networks over non-ideal inter-sensor wireless channels” in Proc. of Int. 
Conf. on Dist. Computing in Sensor Systems and Workshops, Barcelona, 
Spain, June 2011, pp. 1-6. 

[4] R. Abdolee, B. Champagne, and A. H. Sayed, “Diffusion LMS strategies 
for parameter estimation over fading wireless channels,” in Proc. of IEEE 
Int. Conf. on Communications (ICC), Budapest, Hungary, June 2013, pp. 
1926-1930. 

[5] A. Bertrand, M. Moonen, and A. H. Sayed, “Diffusion bias-compensated 
RLS estimation over adaptive networks,” IEEE Trans, on Signal Pro¬ 
cessing, pp. 5212-5224, Nov. 2011. 

[6] R. Abdolee, B. Champagne, and A. Sayed, “Estimation of space-time 
varying parameters using a diffusion LMS algorithm,” IEEE Trans, on 
Signal Processing, vol. 62, no. 2, pp. 403- 418, Jan. 2014. 

[7] X. Zhao and A. H. Sayed, “Performance limits for distributed estimation 
over LMS adaptive networks,” IEEE Trans, on Signal Processing, 
vol. 60, no. 10, pp. 5107-5124, Oct. 2012. 

[8] P. V. Nikitin and K. S. Rao, “Theory and measurement of backscattering 
from RFID tags,” IEEE Antennas and Propagation Magazine, vol. 48, 
no. 6, pp. 212-218, Dec. 2006. 

[9] A. Leon-Garcia, Probability and Random Processes for Electrical En¬ 
gineering. 2nd edition, Addison-Wesley, Reading, MA, 1994. 

[10] A. H. Sayed, Adaptive Filters. Wiley, NJ, 2008. 

[11] R. Abdolee, B. Champagne, and A. H. Sayed, “A diffusion LMS 
strategy for parameter estimation in noisy regressor applications,” in 
Proc. European Signal Processing Conference (EUSIPCO), Bucharest, 
Romania, Aug. 2012, pp. 749-753. 

[12] W. X. Zheng, “A least-squares based algorithm for FIR filtering with 
noisy data,” in Proc. IEEE International Symposium on Circuits and 
Systems, May 2003, pp. IV-444—IV-447. 





























